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Radiation pressure effects in ionization by short linearly-polarized laser pulses are investigated 
in the framework of strong-held approximation, in both nonrelativistic and relativistic formula¬ 
tions. Differences between both approaches are discussed, and retardation and recoil corrections 
are defined. It is demonstrated how these corrections can be incorporated into the nonrelativistic 
approach, leading to the so-called quasi-relativistic formulation. These three approaches are further 
applied to the analysis of signatures of radiation pressure in energy-angular distributions of pho¬ 
toelectrons. It is demonstrated that, for TiiSapphire laser pulses of intensities up to 10^®W/cm^, 
predictions of the quasi-relativistic formulation agree well with those of the full relativistic one, and 
that the recoil corrections contribute predominantly to the radiation pressure effects. 
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I. INTRODUCTION 

Effects in strong-field ionization of atoms attributed 
to the radiation pressure [l| have been studied recently 
both experimentally and theoretically . It is well- 
known that these effects are related to the photon mo¬ 
mentum being transferred to the atomic system. Hence, 
modifying the distributions of ionized electrons. Such a 
phenomenon is absent in the nonrelativistic theoretical 
approach because of the dipole approximation applied to 
the laser field. 

Except for the Ivanov’s paper Q, where the numerical 
solution of the Dirac equation coupled to the 2-cycle laser 
pulse has been considered, in Refs. different types of 
the Strong-Field Approximation (SFA) have been applied 
that go beyond the dipole approximation. Not surpris¬ 
ingly, in all cases the radiation pressure effects have be¬ 
come visible. The aim of this paper is to further develop 
the SFA approach and to analyze angular distributions 
of photoelectrons ionized by linearly-polarized few-cycle 
pulses. This is in contrast to previous works [H where 
circularly-polarized fields were considered in order to de¬ 
scribe multiphoton ionization. Our numerical analysis 
concerns the multiphoton ionization of hydrogen atoms, 
even though we expect that similar effects should be ob¬ 
served for heavier atomic systems. 

This paper is organized as follows. In Sec. m we 
will define properties of linearly-polarized and short laser 
pulses. For such pulses, in Sec. lIIIl we will formulate three 
types of SFA approaches to ionization, each in the ve¬ 
locity gauge. In Secs. [Tv] El and ED we shall compare 
predictions of these approaches and investigate effects re¬ 
lated to the radiation pressure in angular distributions of 
photoelectrons. In Sec. EID we will draw some conclu¬ 
sions and formulate prospects for further investigations. 
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In this paper we keep h = 1. Hence, the fine-structure 
constant equals a = e^/(47r£oc). In our numerical 
analysis we use relativistic units (rel. units) such that 
h = me = c = 1 where me is the electron rest mass. Fur¬ 
thermore, we denote the product of any two four-vectors 
and 6^ as a • 6 = — a}b^ — 

(/i = 0,1,2,3), where the Einstein summation conven¬ 
tion is used. We employ the Feynman notation = 
7 • a = for the contraction with the Dirac matri¬ 

ces 7 ^ and we use a customary notation u = 1 x^ 7 °, where 
is the Hermitian conjugate of a bispinor u. For the 
four-vectors we use both the contravariant (a°, ,a^,a^) 

and the standard (oq, Oa,, Oy, Oz) = (ao,a) notations. Fi¬ 
nally, we use the light-cone variables. For a given space 
direction determined by a unit vector n (which in this pa¬ 
per denotes the direction of the pulse propagation) and 
for an arbitrary four-vector a, we define these variables 
such that all = n ■ a, a~ = — aH, a+ = (a° -I- all)/2, 

and = a — aHn. Thus, the relativistic scalar product 
is a- 6 = a^b~ -\-a~b'^ — and d^cc = dx+dx^d^a;-*-. 


II. LASER PULSE DESCRIPTION 

While in prior investigations of radiation pressure 
effects in multiphoton ionization the circularly- 

polarized laser fields have been considered, in this work 
we consider the laser field of a linear polarization. We as¬ 
sume that the laser pulse lasts for time Tp, which defines 
its fundamental frequency lo = 2tt jTp (being the smallest 
frequency in the Fourier decomposition of the pulse). We 
also assume that the plane-wave-fronted pulse approxi¬ 
mation can be used, which means that the electromag¬ 
netic field is a function of the phase (j) = k-x = ujt — k-x, 
Q < 4> < 27r. We choose the system of coordinates 
such that the laser pulse propagates in the 2 -direction, 
k = (uj/c)n = {ujlc)ez, and its linear polarization vector 
is along the x-axis. This means that the electric field has 
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the form 


5(</.)=fo/5We., (1) 

where fo > 0. The so-called shape function for the elec¬ 
tromagnetic field, fe{4>)} vanishes for </> < 0 and (j) > 2 tt, 
and have a continuous first derivative for all real tp. In 
our numerical analysis, we use the function with the sine- 
squared envelope, 

fsi(p) = ^osin^(^) sin{Nosc(p), (2) 

where iVosc defines the number of cycles in the pulse, and 
the positive real constant Nq is adjusted such that the 
normalization condition: 

^ I dmw? = 2^^osc (3) 

holds. With this normalization, the time-averaged inten¬ 
sity of the laser pulse equals 

where wl = Nosc^^ is the laser carrier frequency, Irei is 
the relativistic unit of intensity, 

/rel = ^ « 2.3 X 10^9 W/cm^ (5) 

ana 

and the dimensionless and relativistically invariant pa¬ 
rameter, 


lefol 

rrieCLU 


= iV„ 


|ego| 

rrieCUJi, ’ 


( 6 ) 


called sometimes the Ritus parameter [3, [1| (see, also the 
review articles iS) , marks the transition between the 
nonrelativistic <C 1) and ultra-relativistic ^ 1) 

regimes of laser-matter interactions. Next, we introduce 
the function 


fW = -[ d(j)'feicp'), (7) 

Jo 

which also vanishes for </> < 0 and cp > 27r, and it defines 
the electromagnetic vector potential, 

A{(p) = Aof{(p)e^ = —f{(p)e^. (8) 

OJ 

In the following, we will use the relativistic notation such 
that A{(p) = Aof{(p)s, where e = (0, e^) and e-e = —1. In 
the nonrelativistic case when the dipole approximation is 
used, meaning that the vector potential and the electric 
field depend only on time, we will assume (p = = ujt. 

In our analysis, we shall use the right-hand-side system 
of Cartesian coordinates defined by the triad of orthonor¬ 
mal vectors (—e^, ej,, e^,). In this system, we introduce 
the polar 0 and azimuthal $ angles (cf. Fig. [TJ. Thus, 



Figure 1. (Color online) Polar 0 and azimuthal angles in 
the Cartesian coordinate system, defined by the triad of unit 
vectors (—e^, e^, ej,). 


an arbitrary vector p is defined by its length |p| and 
spherical angles (0p, <I>p) such that 

P = IpI cos $p -I- By sin $p) sin 0p -|- cos 0p] (9) 
and 

d^p = p'^d\p\d'^flp, d^rip = sin0pd0pd$p. (10) 

This choice of coordinates is motivated by the symmetry 
of the ionization probability distribution, which in the 
nonrelativistic limit and for the linearly-polarized laser 
field in the x-direction is independent of the azimuthal 
angle <I>p, where p is the momentum of photoelectrons. 
This suggests that any deviations from this symmetry 
can be attributed to the radiation pressure. 


III. SFA THEORIES AND APPROXIMATIONS 

The SFA in strong-held ionization has been introduced 
in the seminal Keldysh paper 0 , and during the last 50 
years it has been formulated in different ways and for 
various purposes (see, e.g.. Refs. [iMi)- In the present 
investigations we choose the SFA in the velocity gauge, 
only because in the prior theoretical studies of the ra¬ 
diation pressure effects exactly this gauge (up to 

the Kramers-Henneberger transformation) has been con¬ 
sidered. An extended discussion of these effects in the 
length gauge is going to be presented in due course. 

In the following Sections, we will formulate three ver¬ 
sions of SFA in the velocity gauge. In Sec. IIII Al the Non¬ 
relativistic SFA (NSFA) will be presented, which then is 
followed by the fully Relativistic SFA (RSFA) introduced 
in Sec. IIII Bl It is possible to reformulate the NSFA such 
that it accounts for the interaction of photoelectrons with 
their parent ions (for instance, in terms of the complex¬ 
time trajectories, e.g., [IM3)- For the RSFA it is not an 
easy task, although some attempts have been already un¬ 
dertaken [^, [2l| (see, also the theoretical analysis for the 
Klein-Gordon equation [l^). For this reason, we will for¬ 
mulate in Sec. IIII Cl the Quasi-Relativistic SFA (QRSFA) 
which computationally is not so demanding and, as being 
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very similar to the NSFA, could be further generalized in 
terms of the complex-time trajectories. One of the main 
purposes of this paper is to compare predictions of the 
QRSFA and the RSFA, and to determine for which laser 
pulse intensities the former approach is applicable. 

A general scheme for the SFA is as follows. We split 
the exact time-dependent Hamiltonian, 

H{t)=Ho + Hi{t), (11) 


and q = (g°, q) with = Eq/c. In the following, we shall 
call such four-component objects the ‘four-vectors’ or the 
‘four-momenta’. One has to remember, however, that in 
general they are not the true relativistic four-vectors as 
they do not transform properly under the Lorentz boosts. 

A. NSFA 


into the static, Hq, and the time-dependent, parts. 

The latter accounts for the laser-matter interaction and 
it vanishes asymptotically (i.e., it vanishes in the remote 
past and future). For Hq, we define the bound state \ipo) 
of energy Eq. Then, the exact probability amplitude for 
the transition from the bound state \iPq) to the scattering 
state I'ippit)), labeled by the asymptotic momentum p, is 
equal to 

A{p) = -i y dt(V'p(t)|i?7(t)|V’o)e"‘^“‘. (12) 

\fpp{t)) is the exact solution for the full Hamiltonian H{t) 
which, for t —>• oo, tends to the specific scattering state 
with the incoming spherical waves defined for the Hamil¬ 
tonian Hq. This state is very difficult to determine, there¬ 
fore some approximations are applied. The SFA consists 
in replacing this exact solution by the so-called Volkov 
solution |23l| . In this approximation, the influence of the 
binding potential is neglected, i.e., the probability ampli¬ 
tude in the SFA becomes, 

^sfa(p) = -i J dt{'il^^\t)\Hi{t)\ipQ)e~'^°K (13) 

Here, the Volkov solution fulfills the Schrbdinger 

equation, 

iati4°)(t)) = [i/kin + Hj{t)]\,p(°\t)), (14) 

where i?kin is that part of Hq which remains after neglect¬ 
ing the binding potential. This is equivalent to treating 
the exact solution \ipp(t)) in the lowest Born approxi¬ 
mation with respect to the binding potential. Such an 
approximation is justified if the kinetic energy of pho¬ 
toelectrons is sufficiently large compared to the binding 
energy Eq, or if the action of the laser pulse exceeds the 
influence of the binding potential for photoelectrons. 

In the following discussion, in order to trace similarities 
between the nonrelativistic and relativistic approaches, 
we shall use the relativistic notation in both cases. For 
this purpose, we rewrite Eq. 0 as follows, 

^SFA(p) = -i J ^^^3 J 


Let us now apply the general SFA scheme to the non¬ 
relativistic ionization described in the velocity gauge. In 
this case, the laser-electron interaction takes the form 

-Hi{x) = i—A{k°x°) ■ V -k (17) 

c uIqC zrrieC 


whereas the Volkov solution is defined as 


^pHx) = -^exp 


—ip • X—1 


2meck^ > 


,feA{ (t)')-p 

sCfc° 


(18) 


Here, V is the quantization volume and p = (j)^,p) with 
P° = p^/(2mec). 

For the linearly-polarized laser pulse defined in Sec. m 
the probability amplitude of ionization under the NSFA 
becomes 


^nsfa(p) = i 


■ Mp) 


{rrieCpY ^2 
2mecfc° 


vV J 


d(j) 


rriecp 


exp 


.nieck^ 

-0 ^ 

.p 


£ ■ pf{4') 


fco 


- 4 > + iG'NSFA(</>) 


(19) 


Note that in order to trace similarities with the full rel¬ 
ativistic treatment, in this and in the following formulas 
we keep the factors meC. In Eq. 0, we have introduced 
the laser-dressed energy (divided by c) of the final elec¬ 
tron, 


f=p° + 


rUeCp 

TTlpCk^ 




along with 
Gnsfa(<(>) 

Here, 


= / d<(.' 


meCpL 


irrieCkP 
2meckP 


p{f{4') - if)) 


( 21 ) 


(/) = y^y'V#(0) (22) 


where 

xV;o(9)e-‘«'", 

(15) 

and, similarly, for (/^). Introducing now the Eourier ex¬ 
pansions [ 234 ^. 



(16) 

[/(</))]^ exp [iGNSFA(</>)] = Y2 (23) 

N= — oo 
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for j = 1,2, we arrive at the expression for the probability 
amplitude of ionization such that 


The free-electron Dirac bispinor fulfills the system 
of linear algebraic equations, 


-4nsfa(p) = i-^^^2^NSFA(p), (24) 


- mec)u 


(+) 

pX 


= 0 , 


(33) 


where 


°° -q°-Nk°)/k° _ ^ 

VnsMp) = i{p0 - qO - Nk0)/k0 


^■P ^(1) 

mecp nsfa.tv 


rUeC^l ^(2) ] 

2meCfcO '^SFA.aJ 


^o(p)- 

(25) 


Based on these expressions, we define the energy-angular 
differential probability distribution of ionization, 


d^PNSFA /meCH\'^ me\p\ 


diSpd^Dp 


/meCfxy 

V 27 r / 


27r 


-|25nsfa(p)|^ 


(26) 


and its dimensionless equivalent 


d^Pi 


7^NSFa(p) = 7^NSFa(^^p, ©p, ^p) = , 


or 


7^nsfa(p) = ^ ( 2 n )^ ’ I^nsfa(p)P, 


(27) 

(28) 


which is the probability distribution in the atomic units. 
The latter being presented in our numerical illustrations. 

In the following Sections, we will derive other versions 
of the SFA and, if it does not create misunderstandings, 
we will omit the subscripts labeling these formulations. 


B. RSFA 

In the relativistic theory the procedure is very similar. 
The laser-electron interaction Hamiltonian in the velocity 
gauge has the form, 

Hi{x) = ec7°4(fc • x), (29) 


and satisfies the normalization condition = ^aa', 

with A = ±1 labeling the spin degrees of freedom (2^. 

For our further purposes we define the B-functions, 

■SpA.A.^^) =^pAV^o(g,Ai), (34) 

^^^7.(9) =4tV^o(g,Ai), (35) 

where ipoix, Ai) is the bound state solution of the Dirac 
equation of energy cg° = Eq < nieC^, with the spin po¬ 
larization Ai. The explicit forms of these solutions for the 
gro und state of the hydrogen-like atoms are presented in 
[27j. Since we arrive at the following expression 

for the relativistic probability amplitude, 


AIrsfa(p, A; Ai) = i< 


I nieC^ 

VEo 


d^q 

(^ 




{meCpf ^(0) 


2p ■ n 


B)pixM)r^-x) . (36) 


Now, by analogy to Eq. (EOl), we introduce the laser- 
dressed electron final momentum. 


p = p + 


■ ITLeCp 

p ■ k 


{f)£-P + 


(rWeC^)^ 

2p ■ k 


if) 


(37) 


which fulhlls the conditions (expressed in the light-cone 
coordinates): p~ = p~ and = P^- This allows us to 
write, 

Sp{x) - q - X =(j)~ - q~)x~^ - (p'‘‘ - q^) ■ 

+ (P^ - q^)x~ + G'RSFA(fc°a::“), (38) 


and the corresponding Volkov solution (note that the su¬ 
perscript (-I-) denotes solutions of the Dirac equation for 
electrons, i.e., for particles of positive energies) becomes 
(see, e.g., Ii,[l3l), 



where 



rriecp 
2p ■ k 



-iSp(x) (-t) 

^ “pA > 

(30) 


Sp{x) = p - X — 
and 


f 

JO 


d(j) 


meCpL 


. p ■ k 


f{4')e-P 


(31) 


p = {Ep/c,p), Ep = \/{cpY + {rUeC^Y- (32) 


where 


Grsfa(<(') = [ d(j)' 
Jo 


mecp 

p ■ k 
{rrisCpY 
2p ■ k 


(/(<(') - (/))e-P 


ifi^) - if)) ■ (39) 


Thus, the integrations over x^ and x-^ lead to the conser¬ 
vation conditions p~ = q~ and p^ = . Their solution 

can be put in the form 


q = Q=p+{q^-p^)n. (40) 


After performing the integration d^q in Eq. ((36ll . we end 
up with the expression for the relativistic probability am- 




































5 


plitude of ionization, 


and in the Volkov solution [H, [2^ , 


-4rsfa(p, A; Ai) = i^ 


2 /*27r 

I 


VE, 


d(/)e‘* 


p Jo 


L p 




Following the same ste£s as above, we introduce the 
Fourier decomposition [2J-[2q, 

OO 

[/((?!>)] ^ exp[iGRSFA((/>)] = ‘^RSFA,Ae“'^‘^: (42) 

N——00 


and arrive at 


AIrsfa(p, A; Ai) = i—^^X) rsfa(p, A; Aj), 


with 


°° Q 2 m(p°-q°-Nk°)/k° _ ^ 

25 RSFA(p,A;Ai) = ■^^po_qO_NkO)/k^ 


N— — C )0 


^ [^TX^RSFA,Af-®pA,Ai(Q) 

I IJT-eCp (2) p(0) /^'il (/\/\\ 

^ 2p- k ^RSFA,N^pX,\i (y)J ■ (44) 


(43) 


Hence, the dimensionless probability distribution of ion¬ 
ization can be calculated [cf. the comment below 
Eq. (USD], 


'^RSFa(p) = O'" 


{mec)‘^p^ 


\p\ 


12 


|2^rsfa(p, A; Ai)|" 


A, Ai=i 


(45) 

where the summation over the hnal, A, and averaging 
over the initial, Aj, electron spin degrees of freedom have 
been carried out. This formula is the one that should be 
compared with its nonrelativistic version (1281) . 


C. QRSFA 


Based on the full relativistic treatment, it is now pos¬ 
sible to modify the nonrelativistic theory so it accounts 
for the radiation pressure effects. In doing so, we shall 
follow the proposal put forward by Ehlotzky [2^ (see, 
also Ref. @ and references therein), and discussed in the 
context of the laser-assisted scattering processes [2^. We 
call this approach the quasi-relativistic SFA, as it only 
accounts for linear terms in the nonrelativistic expansion 
with respect to 1/c. For the QRSFA, the modifications 
of the NSFA are done in two places: in the laser-matter 
interaction Hamiltonian, 


C rUeC 2meC 


'^pHx) = ^exp 


1 

VV' 


pk-x 


-ip • X — 1 


d</.'( 


eA{(j)')-p 

rrieck^ — p ■ k 


{(!>') 


2(rneck^ — p ■ k) 


■ k)) 


(47) 


Performing the same analysis as in the last two Sections, 
we can derive the expression for the differential proba¬ 
bility distribution being analogous to Eqs. (1^ and (HSl) . 
The relevant formulas are: 


p = p + 


■ mecp{f)e-p (meC/r)^(/^) 

- rrieck^ — p ■ k 2{mf,ckP — p ■ k) 


k 


(48) 


and 


Gqrsfa(0) = [ d(j)' 
Jo 


■ mecp{f{(j)) - {f))e-p 
meckP — p ■ k 


{meCpY{P{cj)) - (p)) - 
2(mecfc° — p ■ k) 


(49) 


with the corresponding Fourier coefficients jy 

Moreover, we obtain 


OO 

2 ^QRSFa(p) = y]] 

N — — 00 


^ 2 -Ki(p°-q°-Nk°)/k° _ I 

i(pO — Y — Nk^)/kP 


^(2) 

'^QRSFA.AT^ ■ 7* ^eC/rGqj^gp^^^ 


TTleCk^ 

X Vio(p-h (g° -P°)n), 


2mf.ckP 


(50) 


and the dimensionless probability distribution of ioniza¬ 
tion (in atomic units) becomes 

^QRSFa(p) = • |2^QRSFa(p)P- (51) 

We will omit the derivation of the above expressions but, 
instead, we will discuss the physical nature of the substi¬ 
tutions (1461) and (lT7)l . 

The first type of substitutions consists in modifying the 
electromagnetic fields and the vector potential such that 
they depend on both time and space variables. For the 
purpose of the current investigations, we call these the 
retardation corrections. It appears that the only modifi¬ 
cation imposed by these corrections is to shift the argu¬ 
ment of the Fourier transform of the bound-state wave- 
function '00 (p) in (1251) by the vector proportional to the 
unit vector n and independent of the length of photon 
momentum [cf. Eq. (l50l) ]. 

p — >p+{q°-pYn. (52) 


The second type of substitutions, that are present only 
in the electron final state and for the purpose of our fur¬ 
ther discussion we call them the recoil corrections., con¬ 
sists in modifying the denominators in (1471) such that (cf. 
Ref. [13) 

nieck^ —> Wecfc® — p ■ k K, p ■ k. (53) 
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Figure 2. (Color online) Azimuthal distribution of electrons 
ionized by the laser pulse defined in Sec.||T]with Nose. = 4 cy¬ 
cles and wl = A^oseW = 1.55 eV. The time-averaged intensity 
equals I — 2x 10^“^ W/cm^ and the polar angle of the electron 
ejection is 0p = 0.04 757r. The final electron energy is chosen 
to be Ep = rrieC^ + 30 eV for the RSFA and I2me = 30 eV 
for the NSFA and QRSFA. The thick blue circle centered in 
the origin of the coordinate system corresponds to the non- 
relativistic SFA, the thin black curve represents the results 
obtained within the RSFA, and the dashed thin red curve is 
for the QRSFA. 

This is equivalent to introducing into the nonrelativistic 
Keldysh theory the effective momentum-dependent mass, 

me — meff = me - p ■ n/c, (54) 

which accounts for the electron recoil during the absorp¬ 
tion of laser photons. Our numerical investigations pre¬ 
sented below show that this is indeed the dominant con¬ 
tribution to the radiation pressure effects and that, in 
principle, the retardation corrections can be neglected in 
the QRSFA formulation of ionization, as they marginally 
contribute to the overall effect. 

Here, we also note that the recoil corrections are not 
introduced into the interaction Hamiltonian (1461) for two 
reasons: firstly, they would lead to very small modifi¬ 
cations of differential probability distributions and, sec¬ 
ondly, such corrections would not appear in the length 
gauge since both the interaction Hamiltonian —e£{(j)) ■ r 
and the gauge factor exp{ieA{(j)) ■ r) do not depend on 
the electron rest mass. 


IV. AZIMUTHAL DISTRIBUTION 

As we have already mentioned, in the reference frame 
defined in Fig. [T] the ionization probability distribution 
in the NSFA does not depend on the azimuthal angle <l>p. 



Figure 3. (Color online) The same as in Fig. [2] but for I = 
9 X 10^® W/cm^, 0p = 0.027r, and either Ep = rUeC^ -|-200eV 
or 12me = 200 eV. 


This is illustrated in Fig. [51 in which the NSFA is repre¬ 
sented by the thick blue circle with the center located at 
the origin of the coordinate system. The results obtained 
in the RSFA and QRSFA are represented by the nearly 
identical closed curves (almost circles in their shapes) but 
with the center shifted towards the positive z-direction. 
This shift is the signature of the radiation pressure. We 
also observe that for the relativistic and quasi-relativistic 
SFA the results only marginally differ from each other. 
This means that, for the considered intensity of the laser 
field, it is not necessary to apply the full relativistic the¬ 
ory. In this case, the quasi-relativistic SFA very well 
describes the ionization distribution and effects related 
to the non-vanishing photon momentum. Therefore, one 
wonders at which laser intensities the QRSFA fails, show¬ 
ing significant deviations from the RSFA. To answer this 
question, in Fig. Owe compare all three approximations 
for much larger time-averaged intensity and a higher final 
kinetic energy of photoelectrons. As expected, the NSFA 
shows the rotational symmetry but both the QRSFA and 
RSFA exhibit much more complex structures, with many 
lobes that result from quantum interferences induced by 
the radiation pressure. It is also expected that the inter¬ 
ference of probability amplitudes can enhance or suppress 
this effect, in the sense that the radiation pressure can 
either increase or decrease the photoelectron momentum 
along the propagation direction of the pulse. In the two 
figures presented in this Section, we have chosen the polar 
angles 0p (for a given kinetic energy of photoelectrons) 
such that the averaged electron momentum parallel to the 
laser pulse propagation direction, i.e., to the z-direction, 
is positive. This agrees with our intuitive understanding 
of this phenomenon. However, with other choices of the 
polar angle 0p, the radiation pressure effect could be re¬ 
versed and the probability distribution could be shifted 
towards negative z. Such possibilities will be discussed 
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Figure 4. (Color online) The same as in Fig. [3] but omitting 
the NSFA results. In addition, the QRSFA (red dashed line) 
has been modified according to Eq. ([5^. Higher recoil correc¬ 
tions, or even the full relativistic substitution rrieC — >■ p ■ k/k^ 
in the final electron state, do not modify this figure. 


in Sec. Ivll 

A closer look at Figs. [5] and [3] shows that on the y-axis 
the NSFA and QRSFA give almost the same results. If we 
analyze the corresponding expressions for the probability 
distributions we notice that for this particular direction, 
for which n ■ p = 0, they differ only by the momen¬ 
tum shift in the Fourier transform of the bound state 
wavefunction, i.e., by the retardation corrections. This 
suggests that the contributions to the radiation pressure 
effects originate mainly from the recoil corrections dis¬ 
cussed in Sec. IIII Cl In order to account for them prop¬ 
erly it is not sufficient (or even not necessary) to replace 
A{u!t) in the NSFA by A{u!t—k ■ x), but modifications of 
the electron mass induced by the recoil corrections have 
to be considered. 

Further improvements of the QRSFA can be achieved 
by accounting for higher relativistic corrections to the 
effective mass in the electron final state, that go beyond 
the original Nordsieck substitution (1531) . 

rrieC —>■ rriesc = rUeC — p ■ n + - -. (55) 

2meC 

In Fig. S] we compare the full relativistic theory with this 
version of the QRSFA. We observe a very good agreement 
between the RSFA and QRSFA results, which is by far 
better than shown in Fig. [31 Note, that the retardation 
corrections contribute marginally to the QRSFA proba¬ 
bility distribution and, for the considered in this paper 
laser pulse frequencies and intensities, they can be safely 
disregarded. For lower field intensities, the above extra 
modification of the QRSFA can be also neglected. 
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Figure 5. (Color online) Color mapping of the polar- 
azimnthal distribution of ionization, treated in the RSFA, for 
the same parameters as in Fig. O The QRSFA gives the same 
results whereas, for the NSFA, we observe the vertical straight 
stripes (i.e., the distribution is "Fp-independent). The bend¬ 
ing of these stripes is the signature of the radiation pressure. 


V. POLAR-AZIMUTHAL DISTRIBUTION 

A more detailed insight into the radiation pressure ef¬ 
fects can be gained by investigating the polar-azimuthal 
angle distributions of photoelectrons at a given final en¬ 
ergy. When choosing this energy one has to remember 
about the domain of validity of the SFA, as discussed 
above. Note also that, for smaller laser intensities, one 
cannot select too large energies of final electrons since 
the corresponding probability distributions are very small 
and so they contribute marginally to the overall effect. 
For this reason, for the intensity 2 x lO^^W/cm^ we 
choose the electron kinetic energy of 30eV, for which one 
could have doubts about the applicability of the Born 
approximation for the final scattering state of electrons. 
Nevertheless, it can provide an intuitive understanding 
of the radiation pressure effects. In Fig. [5] we present the 
color mapping of the polar-azimuthal distribution for the 
RSFA. Note, that for the QRSFA the map is almost iden¬ 
tical, whereas for the NSFA the bended stripes become 
straight, as for the nonrelativistic case the corresponding 
probability distribution is $p-independent. Because of 
the bended-type character of the relativistic distribution 
it is obvious now that, by choosing properly the polar 
angle of electron emission 0p, the laser pulse can exert 
on the photoelectron either the ‘positive’ pressure, when 
the azimuthal probability distribution is shifted towards 
the positive z-axis (as it has been presented in Fig. 12), 
or the ‘negative’ pressure, when the azimuthal distribu¬ 
tion is shifted in the opposite direction, for instance for 
0p = O.OTtt. This problem is going to be quantified in 
the next Section. 

As for the azimuthal distribution discussed in the pre¬ 
vious Section, for a larger pulse intensity and a higher 
electron kinetic energy, the polar-azimuthal distributions 
become much richer, as illustrated in Figs. [5] and [T] This 
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Figure 6. (Color online) The same as in Fig. [5] but for the 
QRSFA and for parameters dehned in the caption of Fig. |3 




Figure 7. (Color online) The same as in Fig. [6] but for the 
RSFA. 
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time we present the results for both the QRSFA and 
RSFA. Although, in details, these two distributions dif¬ 
fer from each other, they do show astonishing similarities. 
This indicates that, even for the time-averaged intensity 
equal to 9 X 10^^W/cm^, the quasi-relativistic theory suf¬ 
ficiently well describes the radiation pressure effects. As 
before, the more intense laser pulses can exert on elec¬ 
trons both ‘positive’ and ’negative’ pressure depending 
on the selected value of the polar angle 0p. 


VI. AZIMUTHAL-ANGLE-INTEGRATED 
DISTRIBUTION 

The triply differential probability distribution 
7^(i?p, 0p, $p), which has been derived above, can be 
used to determine the averaged momentum of ejected 
electrons, 

(P) = [ 3 d^OppT’(F;p,0p,$p). (56) 

For the chosen polarization and propagation direction of 
the driving pulse, this probability distribution fulfills the 
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Figure 8. (Color online) The probability distribution (I60II of 
measuring the electron momentum along the pulse propaga¬ 
tion direction. The parameters are the same as in Fig. [5l The 
solid blue line corresponds to the RSFA and the dashed red 
line to the QRSFA. 


symmetry condition, 

iP(Ap,0p,27r-$p) =iP(F;p,0p,$p). (57) 

This leads to the conclusion that the y-component of the 
averaged momentum has to vanish, {py) = 0. Since we 
consider a finite pulse, the shape function of the electro¬ 
magnetic potential f{4>) has, in general, a non-vanishing 
constant term in its Fourier decomposition. This means 
that the x-component of the averaged momentum {p^) 
may not vanish. This effect cannot be attributed to the 
radiation pressure, as it also appears in the NSFA. This 
is in contrast to the non-vanishing z-component of the 
averaged momentum (pz)- The latter, in the light-cone 
coordinates denoted as (pz) = (p'l), is equal to 

(p") = / 2 ^^ ^ d^QpP(F;p,0p,$p)(-|p|sin0pCOS$p). 

J OL UIqC 

(58) 

It can be represented as the integral 

(P")= / ^^cl0p|p|tUll(£;p,0p), (59) 

with the two-dimensional distribution, 

FFll(Ap,0p) = -sin^0p / d$p cos$pT’(£’p, 0p, $p). 

Jo 

m 

The advantage of introducing this type of distribution is 
that it is equal to 0 for the NSFA, when no radiation pres¬ 
sure effects are expected. In addition, it quantifies the 
notions of the ‘positive’ and ‘negative’ radiation pressure 
discussed above, as they are strictly related to the sign 
ofIUll(Ap,0p). 

In Fig. [51 we present IVH (Ap, 0p) for Ep = mgCp+30eV 
(for RSFA), p^/2me = 30eV (for QRSFA), and for the 
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Figure 9. (Color online) The same as in Fig.|8]but for param¬ 
eters defined in the caption of Fig.O For larger intensities we 
observe discrepancies between the RSFA and QRSFA, as com¬ 
pared to the results presented in Fig. |5] Note that with the 
extra substitution ISSJ these discrepancies become marginally 
small. 


pulse parameters considered in Fig.[S] As expected, these 
distributions are almost identical. Although, predomi¬ 
nantly, they are positive (i.e., the laser pulse exerts the 
‘positive’ pressure on electrons), they also take negative 
values. If, however, we integrate W^^{Ep,Qp) over the 
polar angle 0p we obtain the positive value. This means 
that electrons of this particular energy are pushed in the 
positive z-direction by the pulse. A similar pattern is 
observed for larger pulse intensities and higher electron 
energies (see, e.g., Fig. ED, but now there are more oscil¬ 
lations of FFII {Ep, 0p) and the discrepancies between the 
RSFA and QRSFA results are clearly visible. Despite the 
obvious differences, qualitatively both distributions look 
similar. Note that, even though for different energies 
and polar angles the laser pulse can affect the electron 
distribution differently (i.e., it can exert either ’positive’ 
or ’negative’ pressure on the electron), its action on the 


whole system ’nucleus-l-electron’ is always ’positive’. 

VII. CONCLUSIONS 

We have studied the effects related to the radiation 
pressure exerted by short, linearly-polarized laser pulses 
on hydrogen atoms during their ionization. In doing so, 
we have presented three types of the SFA developed for 
the velocity gauge. We have shown that, for the Ti- 
Sapphire laser pulses of intensities up to 10^®W/cm^, the 
quasi-relativistic SFA predicts nearly identical effects as 
the full relativistic SFA for the Dirac equation. In par¬ 
ticular, we have discussed modifications related to the 
retardation and recoil corrections. We have shown that 
the dominant contribution to the radiation pressure orig¬ 
inates from the recoil corrections to the nonrelativistic 
SFA. These findings create a possibility to modify the 
complex-time trajectory method for ionization such that 
it can account for the radiation pressure effects. Also, we 
have demonstrated that, depending on photoelectrons fi¬ 
nal energies and polar angles, the radiation pressure ef¬ 
fects can be of ‘positive’ or ‘negative’ types. They are 
more pronounced in the energy-angular distributions of 
ejected electrons than in the integrated, parallel (with 
respect to the laser pulse propagation direction) momen¬ 
tum distributions, where the overall effect is significantly 
suppressed. Extensions of the theoretical methods con¬ 
sidered here to the length gauge and arbitrary laser pulse 
polarizations are under development and are going to be 
presented in due course. 

To summarize, we have demonstrated that, for short 
laser pulses of time-averaged intensities not larger than 
lO^^W/cm^ with the carrier frequency of I.55eV, the ra¬ 
diation pressure effects in ionization can be well described 
by the nonrelativistic Keldysh theory provided that the 
electron mass is replaced by the effective momentum- 
dependent one. The latter accounts for the recoil of elec¬ 
trons during the absorption of laser photons. This has 
been shown by analyzing differential probability distri¬ 
butions in the nonrelativistic and relativistic strong-held 
approximations in the velocity gauge, but we expect that 
the same conclusion is also valid for the length gauge. 
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